Accepted for Publication in Phys. Rev. B 



Commensurability Effects in Large Josephson Junctions 

Leon Balents 

Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4030 

Steven H. Simon 
Department of Physics, Harvard University, Cambridge, MA 02138 



in 
os 

OS 
S3 

a 

(N 
<N 

> 



Os 
O 

Os 

c3 



S3 
O 

o 



Two types of commensurability effects are possible in a large Josephson junction patterned with columnar 
defects. The first occurs for a periodic array of pins when the mean fluxon spacing (tuned by the magnitude of 
the applied in-plane magnetic field) is a rational fraction of the defect spacing. We show that this effect leads, 
under fairly general conditions, to a mapping of the behavior of the Josephson junction near the commensurate 
field values to that of a zero field junction with an effective Josephson penetration depth. The second occurs 
for more general arrangements of pinning sites, when the orientation of the Josephson vortex lattice (tuned by 
the direction of the applied field) nearly matches the orientation of the defects. We investigate this tilt response 
in the limit of a single Josephson vortex. The results are compared, where possible, to recent experiments. As 
an aside from our main analysis, we prove that, contrary to recent claims in the literature, the critical current 
density vanishes in the thermodynamic limit, even in the presence of (non-pathologically distributed) pinning 
disorder. 



I. INTRODUCTION 

The pursuit of large critical current densities in high 
temperature superconductors in external magnetic fields 
has stimulated many investigations into flux pinning. 
Such pinning, caused by naturally occurring and/or ar- 
tificially induced non-superconducting defects, is neces- 
sary to prevent dissipative flux flow in response to an ap- 
plied current. At low temperatures in three dimensions, 
disorder is expected to induce glassy phases, in which 
the current -voltage relation is strongly non-linear, with 
vanishing linear resistivity. These states foe. point and 
line pins have been named the vortex glassa and Bose 
glass,B respectively. Experimentally, continuous normal- 
glass transitions have been observed for both types of 
disorder Ba 

Quantitative comparison with the theory is, however, 
difficult for several reasons. On the theoretical side, nei- 
ther the critical exponents nor the properties of the low 
temperature phase are reliably known.El Experimentally, 
it is hard to prepare large samples with well characterized 
disorder. It is probably this difficulty which has lead to 
a fairly broad range of measured values for these critical 
exponents. 

Josephson junctions (JJs) offer a promising alternative. 
By orienting a magnetic field parallel to the junction in- 
terface, the experimenter can create Josephson vortices 
in the junction, which then acts in many ways like a 1+1 
dimensional version of the bulk vortex state. Because of 
the sophistication and ease of modern fabrication tech- 
niques, samples may be engineered with almost any con- 
figuration of defects in an extremely controlled way. The 
theoretical situation is also much more secure for 1+1 



cal simulation in 1+1 dimensions are also .more tractable 
than their 2+1 dimensional counterparts. 

Patterned JJs are interesting in theijj-pwn right, and 
have been studied by several authors ffl Here, we dis- 
cuss the case of linear (columnar) defects, focusing on 
commensurability effects tuned by changes in the mag- 
netic field. Because of the relative ease of controlling 
the magnetic field versus other system parameters (e.g. 
temperature, disorder), such effects provide a convenient 
probe of the properties of the J J. 

The remainder of the paper is organized as follows. In 
section [01 we review the properties of large uniform JJs 
and the equations describing their behavior. Section |l| 
describes how defects are incorporated into this model. 
In addition, we prove that, other than for particularly 
pathological cases, the introduction of defects does not 
allow a non-zero average critical current density for the 
J J in the thermodynamic limit, contrary to recent claims 
in the literatureJ§BI Commensurability effects which oc- 
cur when varying the magnitude of the applied magnetic 
field in a periodic defect array are explained in section IV . 
Extensive experiments on such systems have been carried 
out in Refs. [To] and [ll], in good agreement with the the- 
oretical results presented here. A discussion of complex 
commensurate states which occur at 1/2-integral fields 
is included in Appendix A. 

A qualitatively different type of commensurability ef- 
fect occurs upon changing the direction of the magnetic 
field. In section we calculate this tilt response_£pr an 
isolated flux line, where the well-known analogyE3 to a 
one dimensional system is particularly useful. In section 
VI, we conclude with a discussion of the experimental 



dimensional glasses, mainly due to | 
renormalization group calculations.1 



re ability to perform 
Large scale numeri- 



checks on this work and its ramifications for future re- 
search. 



1 



II. LARGE JOSEPHSON JUNCTIONS 

A Joscphson junction is an insulating interface between 
two (usually identical) superconductors. For simplicity, 
we will restrict the discussion to the so-called overlap ge- 
ometry, shown in Fig.l. This is a common experimental 
choice, but by no means unique. Other configurations 
behave similarly, but require different boundary condi- 
tions. 

In a J J, all the quantities of interest may be ex- 
pressed in terms of the gauge-invariant phase difference 

7 = 9\ — 8 2 — 2tt/0o dzA z {z) between the two su- 
perconductors, where 0o — hc/2e is the flux quantum.0 
Here Q\ and Qi are the phases of the BCS order param- 
eters in the top and bottom superconducting slab, A z is 
the z-component of the vector potential. The length d 
is the magnetic thickness of the junction, and is related 
to the actual thickness s by d = s + 2A, where A is the 
London penetration depth. The current density in the 
z-direction has the Josephson form 



] z = Jo sin 7. 



(2.1) 



We will regard the tunneling current density jo a s an 
experimental constant. The magnetic field components 
parallel to the junction plane are 



B x — 



0o dj 
2nd dy 
0o 

2nd dx ' 



(2.2) 



Assuming a static situation in the interior of the junc- 
tion, Eqs.2.1 and 2.2 can be combined with the Maxwell 



equation V x B = 47rj/c to yield 
A^V 2 7 = sin 7 , 



(2.3) 



where V 2 indicates the two dimensional Laplacian in the 
interface plane. The length Xj = \/ c0o/ 87r 2 <i7o is known 



as the Josephson penetration depth. Eq.2^3 is known as 
the (time independent) sine-Gordon equation. 

Assuming the current is injected parallel to the-y-axis, 
the boundary conditions augmenting Eq. 2.3 areo 



dx 
dx 

6*7 

dy 

6*7 

dy 



x=0 



x=W 



y=L 



2ndH y 

4>o 

2ndH y 

00 

2ndH x 

00 

2ndH x 



+ 



4ir 2 dl 
c<j) W 

4ir 2 dl 
cfoW 



(2.4) 



Eq.2.3, in conjunction with Eq. 2.4, can have multi- 
ple solutions. In such cases, the equilibrium situation is 
found by minimizing the free energy 



F = e, 



dxdy { -|V 7 - h|" + X-/ (1 - cos 7 ) \ , (2.5) 



where ej = 0g/167r 3 d is the overall energy scale and h = 
2nd(-Hy,H x )/4> . 

We next briefly summarize the properties of a defect 
free junction.Ej In the absence of an applied magnetic 
field, a small or long junction with L < A,/ can carry a 
current / = LWjosinj, with 7 approximately uniform, 
leading to a critical current I c q — LWjo- This can be 
seen by transforming away the boundary conditions ac- 
cording to 7 -> 7+ I(y - L/2) 2 /(2j LWX 2 j). For large 
junctions, where both W, L 3> Xj, the uniform approxi- 
mation fails. In this limit, the phase (and therefore the 
field and current) decays exponentially to zero in the in- 
terior of the junction. This screening is a property of the 
linearized sine-Gordon equation, 



A 2 V 2 7 w 7- 



(2.6) 



In this regime, I c cx A 2 jo- 
in an applied magnetic field, the critical current is fur- 
ther reduced. For strong fields, H 3> 0o/rfAj, screening 
is negligible, and the phase winds approximately linearly 
transverse to the field axis, e.g. 

1-7+ ^ H v x , for H x = 0, H y > o /dAj. (2.7) 



In this limit, the crit ical current, obtained most simply 
by integrating Eq.2.1, takes the Fraunhofer form 



Ic » joLW 



sm(ndHyW/4>o) 



irdH y W/(t) Q 



(2.8) 



For H > 0o/c?A,/, the field penetrates the junction, 
but not uniformly. Instead, it gathers into "Josephson 
vortices" or "solitons" of total flux 0o, across which 7 
changes by 2n. For large fields, these vortices comprise 
only a weak sinusoidal modulation of the field in the junc- 
tio n (a s can be obtained by perturbation theory from 
Eq.p|). 

As H — > H c u, with H c u = 2(f>o/TT 2 d\j, however, the 
solitons sharpen into objects of definite width Aj, be- 
tween which 7 rapidly decays to a multiple of 2n. For 
(H — H c ij) I H c u <C 1, a kind of critical phenomenon, 
known as a commensurate-incommensurate transition 
(CIT), governs the vanishing density of Joscphson vor- 
tices. In this regime, it is legitimate to treat the solitons 
as weakly interacting, since their average separation is 
much larger than Aj. The scaling of the free energy of 



an isolated Josephson vortex can be read off from Eq.2.5 



By using the exact one-soliton solution to Eq.2.3 



7o(a;)=4tan 1 (exp(x/Aj)), 



(2.9) 



and subtracting the free energy for a configuration with 
no vortex, the full expression including the prefactor is 
found as 



2 



F VO rtex — —8(H/H c u — l)—ej. 



(2.10) 



For H > H c u, of course, inter-soliton interactions must 
be included to prevent the vortices from proliferating. 
In addition, we will also allow for small fluctuations in 
the positions and shapes of the vortices, since sufficiently 
small deformations of 7 cost a free energy less than k B T. 
Doing so leads to (neglecting an unimportant constant 
term) the soliton free energy (for N vortices) 



T N = -N(H/H clJ - l)eL + J dy\ £ - 
+eJ2V(\x n (y)-x n ,(y)\/\j) 



dx n (y) 



dy 



(2.11) 



where e = 8ej/Aj, and V(x) ~ 2%exp(— \) is an expo- 
nentially decaying repulsive interaction with magnitude 
and range of order unity (we have approximated the in- 
teractions as local in y_which is adequate for the dilute 
limit considered here).t3 

To unders tand the vanishing of the vortex density at 
the CIT, Eq. 2.11 is used to estimate the free energy as 
a function of i, the average fluxon separation. At finite 
temperature, an entropic contribution must be included, 
due to the wandering of the solitons along the y axis.tZl 
By equipartition, an individual Joscphson vortex wan- 
ders according to {{x{y) — x(0)) 2 ) ~ ksTy/i. Each time 
the fluxon wanders a transverse distance £, it is caged by 
another soliton, and must reverse direction Naively, this 
constraint excludes half of the available configurations of 
the vortex. Taking the logarithm of the total number of 
configurations gives the entropy per vortex, which there- 
fore is reduced by a constant ~ fee In 2. This leads to 
a free energy cost —TAS ~ {k B T) 2 L / {ei 2 ) for a single 
Josephson vortex. Including this with the ene rgetic con- 
tributions from the first and last terms in Eq. ^.ll gives 
the free energy density (per unit area LW) 



ej_ ej_ l/Xj (k B T) 2 A„ 



A 2 / 



(2.12) 



where r = (H/H c u — 1). For k B T -C ej, the minimum 
free energy is determined by the balance between the first 
and second terms, leading to 



£ T=a ~ Aj In(l/r), 



(2.13) 



or, in terms of the internal average magnetic field, (B) ~ 
0o/(Ajdln(l/r)). For sufficiently small r, or at high tem- 
peratures, however, this law breaks down, due to the ef- 
fects of the entropic term. The true asymptotic behavior 
is obtained by balancing the first and last terms to give 



k^T X., 



(2.14) 



The crossover to this behavior is, however, extremely 
close to the CIT at typical experimental temperatures. 



It occurs when r < r co ~ (/c B T/e,/) 2 / ln 3 (ej/fc B T), which 
for, say, a niobium junction with d w 1500A at T — A.2K, 
requires r 10~ 8 . These thermal effects might, how- 
ever, be more observable near the bulk superconductor- 
normal transition (especially in high temperature super- 
conductors) where the Joscphson coupling energy van- 
ishes (within mean field theory) like yl — T/T c . 



III. JUNCTIONS WITH DEFECTS 

The junctions of Ref. |l^ contain defects, which are 
regions in which the local thickness is substantially in- 
creased. For columnar defects parallel to the y axis, 
Eq.2.1 must be modified to jo — > a(x)jo, with a(x) < 1. 



We have assumed that the patterning is uniform in the 
y direction (columnar defects), so that a is independent 
of y. The tunneling current has a strong (exponential) 
dependence when the local thickness variation Sd > £, 
the superconducting coher enc e length. The expressions 
for the magnetic field (Eq.2.2) also become position de- 
pendent through the factor of 1 /d(x) , but. we will ignore 
this weaker dependence in what follows.H 

The modulation in jo(^) results in a modified sine- 
Gordon equation, 



A 2 V 2 7 = a(x) sin 7. 



(3.1) 



The boundary conditions, Eq. 2.4, are unchanged. 

The modified sine-Gordon equation leads naturally to 
def ect -r elated corrections to the dilute fluxon free energy, 
Eq.2.11. To compute them, we first write the modified 



sine-Gordon free energy corresponding to Eq.3.1 



F = e.j J dxdy l^\Vj - h\ 2 + \j 2 a{x) (1 - cos-f)^ . 

(3.2) 

Inserting j(x,y) = J2i^fo(xi(y)), keeping only the lead- 
ing (one-body) term an d ign oring an additive constant, 
gives a correction to Eq j2.1l| of 



■J dyJ2U(x n (y)), 



where 



U(x) = \ J 



dx' a(x + x') 



(3.3) 



(3.4) 



A j cosh 2 (a;'/Aj) 

It has been pointed out by several authorsiJi that an 
appropriate spatially dependent a(x) can lead to an in- 
crease in the critical current of a J J. In section IV we 



will demonstrate this explicitly for the case of periodic 
defects. In the remainder of this section, however, we 
will demonstrate that, contrary to several claimsatl, the 
critical current density of a large (W, L ^> A,/) JJ with 
columnar defects is always zero in the thermodynamic 
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limitS The critical current density also vanishes for the 
long (W ^> Xj, L <C A,/V4'in-line asymmetrical" geom- 
etry considered in Ref. [|J29 The vanishing of the critical 
current density results from the screening of current in- 
herent in the sine-Gordon equation. Even in zero applied 
magnetic field in a pure junction, the critical current den- 
sity decreases with the system size once it is larger than 
Aj. 

To prove that the critical current density indeed van- 
ishes, we will take a more concrete model of the disorder. 
Specifically, a(x) is taken to be piecewise constant on 
intervals of variable length. On all intervals, we require 
< a < 1, the upper bound corresponding to a perfect 
junction (in fact, it is only necessary that a is bounded 
above and below by arbitrary constants). The lengths 
of the intervals should also be well behaved, so that the 
mean length J2i Xi/N = x exists. For simplicity, we con- 
sider the one dimensional junction (i.e. L <C Aj) studied 
in Ref. | 

Consider the mean current density for the junction, 



1 



j = j dxj a(x)sinj, 



(3.5) 



where 7 is the solution of the one-dimensional sine- 
Gordon equation, 



X J^ = a(x)amj. 



(3.6) 



Multiplying by jo and integrating over x gives 

j=j \ 2 j(v(W)-v(0))/W, (3.7) 

where v[x) = d r y(x)/dx. Multiplying Eq. |3.6| instead by 
Xj 2 v(x) and integrating from x to x' , one finds 



E(x') - E(x) — X 2 ^2 §aiCosj(xi), 



(3.8) 



leads to a similar bound. We have therefore established 
in general that 



IjI < Cj - 



A., 



(3.10) 



which vanishes in the thermodynamic (W — > 00) limit. 

This result can be understood in terms of an amusing 
analogy. Consider a simple pendulum in a time varying 
gravitational field. If the gravitational constant can not 
become negative or exceed it's normal value, Eq.3.6 de- 
scribes this situation with 7 giving the angle of t he pen- 
dulum and x taking the role of time. By Eq. 3.7, a non- 



zero critical current density would be equivalent to saying 
that one could continuously accelerate the pendulum by 
varying the gravitational field. Clearly, as the pendulum 
begins to move faster, the field must be switched from 
strong to weak more and more rapidly to continue to give 
it an average acceleration. This would require a vanish- 
ing x, corresponding to the infinitesimal time interval in 
the large W limit. 



IV. PERIODIC DEFECTS 

Although a finite critical current density is impossible 
for an infinite system, an increased critical current can 
be obtained in finite systems by introducing pinning. A 
natural choice for such pinning is to choose a defect pe- 
riodicity commensurate with the Josephson vortex spac- 
ing in the junction. Such junctions have been studied 
extensively in Refs. |lO| , |ll| , and we discuss a theoretical 
approach to these junctions here. 

For ease of presentation, we will assume that the defect 
configuration is such that the unit cell can be chosen 
even. The extension to unit cells without symmetry is 
str aigh tforward. In this case, the phase is governed by 
Eq.3.1, where a(x) can be written in the form 



where xt is the coordinate of the boundary between the 
i th and (i + l) th constant region, and Scti is the jump 
in a at that boundary. The function E(x) = v(x) 2 /2 + 
a(x) cos 7(2;). Using |<5ai| < 1 and the existence of x, the 
magnitude of the right hand side of Eq.|3.8| is bounded 



above by Aj N = A ; 



\x — x 



lx (for large L). Since the 



cos 7 term in E(x) is order one, we have 
\v{x') 2 ~v{xf\ < CXj 2 \x-x'\/a 



(3.9) 



Here and in the remainder of this section, C indicates any 
dimensionless constant of order one. Consider first the 
case where v(W) and v(0) have the same sign. Then 
\v(W) 2 - v(0) 2 \ = \(v(W) - v(0)){v(W) + v(0))\ > 
\v(W) — v(0)\ 2 . Using Eq.3/7, this gives the bound 
< CjoXj/^/xW. If v(W) and v(0) have opposite 
signs, then there exists some x* between and W with 
v(x*) = 0. Both x* < W and W — x* < W, so 
\v{W) - v(0)\ = \v(W)\ + \v(0)\ < 2C\j 2 W/x, which 



a [x) 



E 

m=0 



a m cos mqx, 



(4.1) 



where q = 2ir/a, with defect lattice spacing a. Ma king 
the change of variables 7 = nqx + r] and inserting Eq.4.1 
Eq.|3.1| becomes^ 



X 2 ^ 



m— 



sin [{n + m)qx + rj\ 



+ sin [(n — m)qx + rj] 



(4.2) 



Under this tran sformation, the form of the boundary con- 
ditions (Eq.2T) remains unchanged (with rj replacing 7), 
except that the y component of the magnetic field is re- 
placed by an effective value 



H v - H y - 



nq4>0 
2nd ' 



(4.3) 
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To simplify further, we will assume that the defect lattice 
is not too large, so a < Aj, and that the field has been 
chosen close to commensurate, so that H y d\j -C </)$. Un- 
der these assumptions, r\ will be slowly varying, and all 
the sines but the second m = n term will average outEHl 
(u nless n — 0, in which case both m = terms are equal). 
Eq. 4.2 then reduces to a uniform sine-Gordon equation, 



sin 77, 



(4.4) 



where the effective Josephson penetration depth is 



I y/l/a Xj 



n ^ 
n = 



(4.5) 



Although we have not discussed time dependence in any 
detail, exactly the same mapping holds if it is included 
through th e us ual X 2 j(cd 2 + /39t)7 terms on the left hand 
side of Eq.2.3. Thus the statics and dynamics near the 
commensurate fields will behave as in a uniform junc- 
tion, with an n-dependent effective penetration depth 
and shifted magnetic field. For small H y , the physics 
will therefore extremely closely mimic that of the uni- 
form case. In particular, the critical current at the n th 
peak is simply evaluated at zero field. Thus 



f(L/\j)LWj 



f(L/\j)LWj ^, 



(4.6) 



where /(%) -> 1 for % < 1 and f(x) ~ 2/x for x > 1- 

As an example, we explicitly compute the effective 
parameters for the case studied experimentally in Ref. 
|lO| . There, the thickness of the defects is large enough 
(d 3> £) to render a — on the pinning sites. In this 
case, the function a(x) is a periodic train of rectangular 
pulses of width a — Wd and height 1, where Wd is the de- 
fect width, separated by spaces of width Wd with a = 0. 
Choosing th e or igin at the center of a defect, the Fourier 
series in Eq.4.1 can be inverted to yield 



_ ( -2sm(imwd/a)/(irn), n>0 , ? , 

an -\l-Wd/a, n = ■ 

This implies an effective Josephson penetration depth of 



A, = 



J Xjy/trn/ sm.{imwd/ a), 



Aj/V 1 - i»d/a, 



n> 
n = 



(4.8) 



For an effectively long or small junction (note that this 
only requires Aj > L, not Xj > L), this gives the crit- 
ical current I c = LWjo\ sm(Tmwd/a)/(irn)\ for n > 
and I c — LWjo(l — Wd/a) for n = 0. In this limit, 
these results have been obtained in Ref. [To] using a sim- 
ple current blocking model. We point out here that 
the region of validity of these results is actually much 
wider than would be naively predicted, because the in- 
creased effective Josephson penetration depth brings the 
system further into the small junction limit. The pre- 
dicted modification of the commensurate critical current 



due to th e in crea sed effective screening length (obtained 
from Eqs.4J3 and- 4^) has recently received strong exper- 

tJST 



imental support! 

Further, we can discuss in this way the width of the 
commensurate peaks. For Aj <C W , in the overlap geom- 
etry considered here, the width of a peak is determined 
by the screening condition H y ~ H c \j ~ <po/dXj. Using 
Eq.4.8, this gives 



AH: 



screened 



$j y/sm(Tmw d /a)/mr, 
^-^Jl-Wd/a, 



n>0 
n = 



(4.9) 



We note in passing that an in-line geometry, in which 
the current is injected transverse to the magnetic field, 
has considerably different behavior in this limit, due to 
the restriction of the current to screening layers near x = 
0, W . In the unscreened limit (still the overlap geometry) 
where Xj » W, 



ah: 



unscreened 



go 

dW 



(4.10) 



Between the two scaling regimes there is presumably a 
smooth crossover (which could, in principle, be calculated 
from a solution of the uniform sine-Gordon equation). 



V. TILT RESPONSE OF A SINGLE FLUX LINE 

As discussed in section ||, the fiuxon density vanishes 
continuously as H is reduced toward H c \j. Once the 
mean spacing I > 10Aj, say, the exponentially decay- 
ing interactions between Josephson vortices are weak, 
and it is reasonable to attempt to approximate the soli- 
tons as independent. In fact, at finite temperature in 
an infinitely long (L — 00) system, thermally excited 
transverse wandering of the vortices actually invalidates 
this approximation. In normal experimental conditions 
away from the bulk superconductor-normal transition, 
ho wever , k^T/ej is so small (c.f. the paragraph after 
Eq.2.14) that this fluctuation effect is negligible for typ- 
ical sample lengths. A more serious problem is the ex- 
tremely weak divergence of I near H c u given in Eq.2.13 



means that it is necessary to reach a reduced magnetic 



field of r < e 



5 x 10 to get to this regime. 



It is nevertheless interesting from a theoretical point of 
view to consider the limit of an isolated flux line. Such 
a situation is somewhat more realizable in three dimen- 
sions in a bulk superconductor, a case into which we may 
hope to gain insight through this simpler model. It also 
represents a soluble limit, even in the case of a random 
distribution of defects. We discuss this interesting ran- 
dom case next at T = 0. 

For simplicity, we define the slop e of t he appli ed m ag- 
netic field h 



Combining Eq. 2.11 and Eq.3.3, the 



free energy for a single flux line is given by 



5 



T = 



dy 



dx(y) 



dy 



- h 



U(x(y)) 



(5.1) 



-h(x(L)-x(0)) + / dy 
Jo 



dx{y) 



dy 



U(x(y)) 



constant 



(5.2) 



where we have set the overall energy scale e = 1 for 
simplicity. Since we are neglecting thermal effects, the 
physics is completely independent of this choice. 

This free energy is also the expression for the classi- 
cal action of a particle of mass 1 in 1+1 dimensions in a 
potential — U[x] where y corresponds to the time coordi- 
nate for the classical particle. Note that this mapping is 
somewhat different from the usual boson analogy.E£l Since 
—h(x(L) — x(0)) is a boundary term, it does not effect 
the equations of motion for the classical particle. The 
quantity 



E=~ 
2 



dx 
dy 



- U{x{y)). 



(5.3) 



corresponds to the energy of the classical particle system 
and is thus conserved. This means that the local slope of 
the flux line ^ is a function only of the local potential 

U[x] and of a single constant of the motion E. 

The density of states for the flux line potential U is 



D(u) 



I rW/2 
lim — - I dx 8(u 

V^oo W J_ W/2 



U\x] 



(5.4) 



For now, we will require that the potential U is bounded 
from below, which means there exists a 



min U > — co. 



(5.5) 



Note that D(u) =0 for u < u min . 

We will compute the tilt response, defined as the slope 



g ^ _ x(L;h)-x(0;h) 



(5.6) 



of the flux line in its ground state. For small angles, 
9 w (B x ) I {By). If h = 0, the minimum energy solution is 
x(y) = xq where C/[xo] = Umin, i-e. the flux line in pinned 
to the deepest point of the random potential. This yields 
a free energy of 

F = u minL (5-7) 

If the magnetic field is applied at a nonzero slope h, it is 
possible that the flux line remains pinned on the deepest 
point, yielding the exact same value for the Hamiltonian. 
However, it is also possible that the flux line becomes 
unpinned and traverses the system at some average slope 
9 defined as above. Using Eq.5.3, the magnitude of the 
local slope is given by 



dx 
dy 



E + U(x(y)) 



(5. 



It is a simple exercise to show that, if the upper end- 
point x(L) is free and the lower endpoint is at an ab- 
solute minimum of U(x), the minimal path with always 
have dx/dy > for h > 0. It should be noted that, for 
an unpinned path (dx/dy > 0), E > —u min , so that the 
argument of the square root is nonnegative. In the lan- 
guage of the classical particle this is just the requirement 
that the total energy of the particle is larger than the 
potential of the highest hill, so it has e nough energy to 



traverse the system. Integrating Eq.5.8. the inverse slope 
o 1 ' 



8 x(L)-x(0) 



IS 



a(E) 



1 



" (£) . dy 
dx — 

x{L) - x(0) J x(0) dx' 



(5.9) 



In the limit of a large system when the slope 9 is nonzero, 
the flux line samples all of the disorder in the potential 
equally. Thus since the probability that a given point x 
has potential U[x] = u is just the density of states D(u), 
Eqs.5^ and give 



ct(E) = du 



D(u) 



2[E + u 



(5.10) 



so long as 9 ^ 0. This relation implicitly determines E as 
a function a. The slope is determined by optimizing over 
all possible values of E. It should be noted that Eq. 5.1C 
is completely independent of any spatial correlations in 
the function U(x) so long as the system is large enough 
that the flux line samples all of the disorder equally. 

By using Eq.|]^ to eliminate U(x), Eq.|j can be 
rewritten, up to an additive constant, as 

f x{ - L) dx 

T= -h(x(L) -x(0)) -EL+ / dx— , (5.11) 

Jx(o) d v 



Taking L — > co, 



F{E,h) = L 



a(E) 



E 



a(E) 



duD{u) 4/2 I E 



(5.12) 



Extremizing T with respect to E gives 



duD(u) ^2 (e + u^ (5.13) 



which in principle, can be inve rted to yield E(h) and the 
tilt response 9(h) (cf. Eq.j5T(i|). For future use, we note 
in passing that 



dh_ 

dE 



a(E). 



(5.14) 



G 



For s ome values of h, there will be no E satisfying 
Eqf .13[ In this case, there is no finite optimal value 
of the inverse slope a, and the ground state of the flux 
line is pinned with x{y) = xq, and U[xq] = u m in so that 
9 = 1/a = 0. As ft, is increased from zero, the slope 9 
remains zero (pinned) until h reaches some critical value 
h c , where 9 begins to increase (unpinned) to nonzero val- 
ues. We therefore have dB x /dH x \} l= o — 0, i.e. there is a 
transverse Meissner effect. 

The behavior near h c constitutes a critical phe- 
nomenon analogous to that at the low er cr itical field . 
Using the optimal value of E found in Eq. 5.13 in Eq. [5.12 
the optimal free energy for a given h is 



mm^iE, h) = -E(h)L 



(5.15) 



as L — ► oo. This free energy for th e flu x line is equal to 
the energy of a pinned flux line (Eq.5.7) when E obtains 
its minimal value —u m ; m . In terms of the language of 
the classical particle, this simply says that the classical 
particle becomes "unpinned" as soon as it has sufficient 
energy E to get over the highest hills of the potential. We 
will be able to expand the necessary equations around the 
critical value E = —u m i n and more closely examine the 
transitional behavior. 



1. An Explicit Example 

To illustrate these results, consider a disorder potential 
with the uniform density of states given by 



D(u) 



-1/Ur, 







otherwise 



(5.16) 



with u m in < 0. Eq.5.1C can be explicitly integrated to 
yield 



V2 



E- \ E + U r , 



(5.17) 



E c = —Umin, so the critical value of the slope of the flux 



As discussed above, the critical value of E is given by 
E c = - 
line is 

1 



1 



a{E c ) 



(5.18) 



Similarly, Eq.5.13 can be explicitly integrated to give 



HE) 



-3u r , 



E 3 ' 2 - (E + um m f' 2 (5.19) 



which yields a critical slope of the applied field 



h c — h(E c ) — -V—2umin- 



(5.20) 



For slopes of the applied field that are less than h c , the 
flux line is pinned at 9 = 0. When h reaches h c , 9 jumps 



to 9 C , and then increases continuousl y as h is further 
increased. In this particular case, Eq. 5.17 can also be 
explicitly inverted to yield 



E(a) = ~ 



V2 



-1 2 



V2 



(5.21) 



By substituting Eq.5.21 into Eq.5.19, the complete (al- 
beit messy) analytic form for the function h(a) can be 
found. In general, however, a(E) is not explicitly invert- 
ible, so an analytic form for h(a) can not be obtained. 
Furthermore, the complete function h(pt) is clearly very 
sensitive to the precise form of the density of states. 
Thus, we would like to focus on physics that is in some 
sense more universal. 



2. Angular Exponent Near The Depinning Transition 

To this end we consider the behavior near the critical 
slope. For h — h c small but positive, we should have the 
relation 



(h - h c ) ~{9- 9 C )~ 



(5.22) 



for some exponent 7. As long as 9 C is nonzero, this then 
implies 



(h - h c ) ~ (a - a c y. 



(5.23) 



In general, the value of the exponent 7 will only depend 
on behavior of the density of states D(u) very near to 
Umin- In the above example, using the analytic form for 
h(a), it is trivial to expand around a c to get 7 = 2. In 
fact, 7 will take the same value for all systems where the 
density of states D(u) has a finite jump discontinuity at 

Umin • 

Let us suppose that the density of states takes the form 

D{Umin) ~ (H - Umin) 13 (5-24) 

for u close to but greater than u m in for some exponent 
j3. Of course we must have > —1 such that the integral 
of D is finite and D is normalizable. 

We are now concerned with the behavior for E close 
to but greater than —u m in- We define 

AE = E + u min - (5.25) 

With AE small, Eq.5.10 implies (with p = u — u m in) 



a(E) - / dpp p (p + AE) 
Jo 



-1/2 



(5.26) 



where q is some small cuto ff beyond which D(u) no longer 
has the scaling form 5.24. This integral is finite for (3 > 



-k, and otherwise diverges as 

a(E) - (A£)^ + J' for < -\. 



(5.27) 
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On the other hand, Eq.5.13 implies that h(AE) is always 
finite. But from Eq.5.14, we have dh/dAE ~ a finite as 
AS — > only for (3 > — \ as described above. Finally, 



dot 
dE 



dpp p (p + AE) 



-3/2 



which is finite for (3 > i and diverges as 



^ ~ (AE)V- 
dE 



for /3 < i 



(5.28) 



(5.29) 



In the case of /3 > ^, we have h, a = dh/dE, and 
da/dE all finite and nonzero as AE — > 0. Thus, expand- 
ing around AS = 0, we have (h — h c ) ~ (a — a c ) ~ 
(0- C ) - AS, so 7 = 1, with 7i c > and 9 C = l/a c > 0. 

In the case of h > (3 > — h, dh/dE = a is again finite 

as AE — •> 0, so (h — h c ) ~ AS, but da/dE diverges as 
A _B(/3-i/2) j SQ ( a _ ftc ) „ A £(/3+i/2). Thus 7 = 1/(^+1) 

with /i c > and 8 C > 0. Here we recover the result 7 = 2 
for (3 = that we derived in the explicit example above. 

In the case of— 1 < [3 < — i, a = dh/dE is divergent 
as A^t/3+1/2) as A £ ^ 0j so (ft _ „ A e(P+3/2)^ 



P+3/2 
'/3+1/2' 



and 



Here a c is divergent, so C = 0. Thus, 7 
Eq. |5.23| does not hold. 

For j3 approaching — A the exponent 7 diverges. For 
(3 = — i, Eq. |5.28 can be integrated for small AS to 
yield a ~ — ln(AS) + ■ • • which is divergent for small 
AS so C — 0. Since dh/dE — a, we have (h — /i c ) ~ 
AS(lnAS- 1). Thus, 



(h - h c ) ~ e 



-i/e i 



for (3 



(5.30) 



3. Finite Sized Systems 

In the case of a finite sized system, a flux line with 
nonzero slope can only sample all of the disorder equally 
if \x(L) — x(Q)\ is much larger than the correlation length 
of the disorder. If \x{L) — x(0)| is not much larger than 
the correlation length, we can not use the above method 
to average over the disorder, and we must treat the cor- 
relations explicitly. 

For example, in any finite sized system with smooth 
disorder (ie, U(x) has nonzero correlation length), there 
is some genericly unique point xq such that U(xq) = 
Umin- Near this point, we can expand U(x) as 



U(x) 



+ A(x~x ) 2 + 



(5.31) 



For h = 0, the lowest energy solution is the flux line 
pinned at xo. For nonzero h, the flux line should pull 
away from this minimum in a continuous manner. In 
the simple case of U(x) — Ax 2 , (choosing u m i n = 0, 



and Xo = 0), we can study this behavior explicitly. For 
simplicity we fix one end of the flux line to be pinned in 
the middle of this deepest well, that is, x(0) — xq = 0, 
and we allow x(L) to vary. Then L, S, and x{L) are all 
related by the restriction 



dy [*V 
L = I dx— = I dx 

dx 



x(0) 



1 



x(0) 



2(E+U(x)) 



(5.32) 



which can be integrated explicitly in this parabolic case 
and solved to yield 



S = A[x{L)f 



sinh 



(LV2A) 



(5.33) 



Now differenti atin g Eq. 5.11 with respect to x{L) (and 



again using Eq.5.8 to find a minimum of the f ree energy, 
yields the general result (compare to Eq. 5.13|) 



h = 



E + U{x(L)) 



(5.34) 



which in the case of the parabolic potential yields 



z{L) = h 



/ 2Acoth ( LV2A 



(5.35) 



This calculation can equally well be performed in the case 
where we allow both endpoints to vary and the result is 
similar. Of course this result holds only in the range 
where x(L) is small enough such that the potential U 
still looks parabolic. Furthermore, if x(L) becomes suffi- 
ciently large that there is a point x\ far away such that 
U(xi) < U (x(L)) then it is possible that the flux line will 
stretch to a more favorable position near x\ . Nonetheless, 
for a generic finite sized system, for h sufficiently small, 
the slope 8 = [x(L) ~ x(0)]/L will vary as 8 ~ h. How- 
ever, as can be seen from Eq.5.35, as L gets larger, the 



range of h values for which this parabolic approximation 
is valid is reduced exponentially. For larger values of h or 
for larger L, the system may then cross over to behavior 
more like the infinite sized system discussed above. 



VI. CONCLUSION 

In this paper, we have discussed two types of com- 
mensurability effects in Josephson junctions containing 
columnar defects. The first kind is tuned by changing 
the magnitude of the applied magnetic field along the di- 
rection of the defects. For periodic arrangements of the 
pins, the vortex lattice induced by an applied magnetic 
field locks in to reciprocal lattice vectors which are inte- 
gral multiples of the defect array wavevector (fractionally 
commensurate states also occur, but are suppressed by 
several powers of a/Aj, as discussed in appendix A). Near 
these states, the critical current of a typical junction is 
substantially enhanced, though the critical current den- 
sity remains zero in the thermodynamic limit. 
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A second type of commensurability can be investigated 
by varying the transverse component of the field, or its 
direction. This effect is present for many distributions 
of the defects, though the character of the transition to 
the incommensurate (unaligned) state depends upon this 
distribution. In the case of periodic defects, the map- 
ping of section IV shows that a commensurate vortex 



lattice remains locked in until a finite threshold trans- 
verse field is reached, after which kinked fluxons enable 
the array to tilt. The transverse constitutive relation 
Bj_(H±) was found, under rather general conditions, to 
be simply a rescaled version of the zero parallel field re- 
sult. We also investigated the dilute limit (near H c \j) 
where each vortex acts independently. Using a mapping 
to one-dimensional particle mechanics, we showed that 
in most random distributions of pins, the fluxon makes 
a first order jump to a non-zero tipping angle as the 
transverse field is increased to a critical value. Above 
this field, this angle increases continuously, and near the 
critical field this increase can be described by an almost- 
universal scaling exponent. Similar discussions of this 
"tilt response" have been given for the case of a finite 
density of fluxons with randomly placed defects in Ref. 
|23| . In a periodic array of defects away from commensu- 
ration, and in the single vortex limit, the effect does not 
exist at finite temperatures, as discussed in Ref. ^J. 

Recent experimental work on long overlap junctions 
has confirmed many of our theoretical, conclusions for the 
case of the periodic defect array.E£lEJ In particular, the 
critical currents at the commensurate peaks agree quite 
well with the estimate following Eq j4.§| , with small devi- 
ations which are still well understood within the model. 
Although the peak widths are less reliably obta ined, they 
appear to be consistent with Eqs.[h9 and 4.10. 



A particularly dramatic consequence was reported in 
Ref. [n]. There, the I-V characteristics were measured 
above I c at various commensurate magnetic fields. These 
I-V curves showed "commensurate field steps," analo- 
gous to the zero field steps seen at zero field. Upon 
rescaling the current axis by [Xj(n — 1)/Xj(n = 0)] 2 , 
the steps for n = and n = 1 collapsed to a single 
curve. This dynamic measurement confirms the general- 
ity of the mapping. Another potential characterization is 
the number of zero field steps (which is roughly L/2Xj), 
but an accurate determination requires an impractically 
large system width. It would be interesting to further ex- 
plore experimentally and theoretically the consequences 
of this mapping. 

Another interesting possibility from the theoretical 
point of view is a "floating" phase, in which a commen- 
surate fluxon array is prevented from locking into the 
periodic defect array by thermal fluctuations. As dis- 
cussed in [n[ the large value of ej/(k^,T) renders thermal 
fluctuations unimportant unless 1 — T/T c <C 1. Such a 
floating phase may, however, be possible near the bulk 
superconducting transition, along with other interesting 
fluctuation effects not discussed here. Experimental ob- 
servation of such a state would no doubt be difficult due 



precisely to these noisy fluctuations. We must leave these 
experimental and theoretical questions open. 
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APPENDIX A: HALF-INTEGER 
COMMENSURATE STATES 

In this appendix, we extend the treatment of section 
IV to handle the behavior near the half-integer commen- 
surate fields. This condition is defined (at the center of 
a peak) by 



R (k/2) 

y 



2 2nd' 



(Al) 



where k is an odd integer (for k even this condi- 
tion reduces to that for the integral peaks discussed 
earlier). Higher denominator fractional commensurate 
states should also occur, but the treatment becomes con- 
siderably more complicated and will not be described 
here. 

Proceeding as before, the change of variables 7 = px + 
r], with p = kq/2, leads to the equation 



A 2 V 2 r, 



E 

m— 



sin[(p + mq)x + rj\ 



+ sin[(p — mq)x + rj\ 



(A2) 



Unlike Eq. 4.2, Eq.A2 does not contain any non- 



oscillatory terms. More precisely, a first order pertur- 
bative expansion in a will give convergent results. Nev- 
ertheless, we expect that the fluxon array should lock 
in to the defect lattice for small H v . In fact, the half- 

;e in 
i To 



integer commensurate effect shows up as a div erge 



the second order perturbative solution of Eq.A2 



understand the nature of this divergence, we make the 
(non-linear) change of variables from r\ to a, where 



V = Vo(<y,x) + cr, 



(A3) 



where 



Vo = 



EOl m 
Wr 



sin[(p + mq)x + cr] 



+ 



[p + mq) 2 

sin[(p — mq)x + a] 
(P - mq) 2 



(A4) 







To motivate Eq. A3 , note that 770 is the lowest order solu- 
tion which would have been obtained from Eq. A2 , taking 
a constant value of 77 = a on the right hand side. Despite 
the non-linearity, 



V 2 ?/ « V 2 cr + ' 



sin[(p + mq)x + a] 



+ sin[(p — mq)x + a] 



(A5) 



up to oscillatory terms which may be neglected to lowest 
non-trivial order in 1/JXjp) for small H y . At this level 
of approximation, Eq.A2 becomes 



JJ, we therefore expect only minor quantitative changes 
in behavior. Physically, one effect of these corrections is 
to account for overlap with the long "Fraunhofer" tails 
of the integer-field (including zero-field) commensurate 
peaks. A crude way of accounting for this ef fect is to 
treat the critical current, calculated as in as a 

shift relative to the (experimentally determined) Fraun- 
hofer background. 

Neglecting these complications, and taking into ac- 
count the final change of variables from a to a, the effec- 
tive current and effective fields are 



s'm[(p + mq)x + 7/0 + cr] 



H y = 2{H y -H^) 

H x = 2H X , 
1 = 21, 



(All) 



+ sin[(p - mq)x + 7/0 + a] - sin[(p + mq)x + a] where H v k/ 2) is t he commensurate field defined in Eq.|Al 



sin[(p — mq)x + a] 



(A6) 



From Eq.AlO 



The leading (linear) contribution in a Taylor series in i] Q 
gives the result 



A 2 V 2 cr = £ 



-170 



cos[(p + mq)x + a] 



it is clear that the half-integer states 
have critical currents suppressed by a factor of (a/Aj) 2 . 
We expect that higher denominator fractional states can 
also occur, but are suppressed by even higher powers of 
a/Xj. 

We have estimated the magnitude of the critical cur- 
rent peak in the experiments of Ref. [HI Using Eqs.4.7 



m=0 

+ cos[(p — mq)x + a] 



(A7) 



By inserting Eq.A4 into Eq.A7, combining sines and 
cosines, and neglecting oscillatory terms independent of 
cr, one finds 



A 2 V 2 cr = ]T 



1 



, 8A 2 I (p - m'q) 2 



([(m + m')q — 2p]x 



A10, 4.6, and All, we obtain a lowest 1/2 integral peak 
reduced by a factor of approximately 600 from the zero- 
field critical current. Using I c o = 30m A as determined 
experimentally, this gives AI C w 0.05mA A feature 
is barely visible in the critical current versus field plot 
(Fig.2 of Ref. with a deviation from the background 
of roughly 0.2mA Given the crude treatment of the 
boundary conditions and the smallness of the effect, such 
order of magnitude agreement is satisfactory. 

([(m — m')q + 2p]x + 2cr^ 



2cr 



sin I 



1 



(p + m'qY 



sin I [(m + m')q + 2p]x + 2a) — sin I [(m — m')q — 2p\ — 2a 



(A8) 



Keeping only the unmodulated contributions propor- 
tional to sin(2cr), we arrive at the effective equation 



AjV 2 c> = sin <r, 



(A9) 



where a = 2cr and the effective Josephson penetration 
depth in this case is 



Aj = Aj \4A 2 ^ 



1 



\ ^ 



jP 



4A J n 

J m— 



^mCK|m— k\ 



\p — (k — m)q\ 2 



\p — (k + m)q\ 2 



-1/2 



(A10) 



It should be noted that the change of variables (Eq.jAj) 
from 77 to cr induces a com plic ated transformation of the 
boundary conditions in Eq. [2.4| . The corrections are, how- 
ever, small by a factor of \/{\jp) 2 . Although we have 
not performed a detailed analysis of their effects upon the 
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FIG 1: Overlap geometry for a large Josephson junction. 
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